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Abstract 

Taking proton as an ensemble of quark-gluon Fock states and using the 
principle of detailed balance, we construct a simple statistical model for parton 
distribution of proton. The recent observed Bjorken-x dependent light flavor 



sea quark asymmetry d(x) — u(x) can be well reproduced by Monte Carlo 
. simulation as a pure statistical effect. 



PACS: 12.40.Ee, 12.38.Lg, 14.20.Dh, 14.65.Bt 

Proton is the simplest system in which the three colors of QCD neutralize into a 
colorless object, but its internal quark-gluon structure is still not well understood. 
The complication comes from the presence of sea quarks in the proton. In all global 
analyses of parton distribution in nucleons before 1990, a symmetric light-quark (u, 
d) sea was assumed, based on the usual assumption that the sea of quark- ant iquark 
pairs is produced perturbatively from gluon splitting |T|. However, a surprisingly 
large asymmetry between the u and d sea quark distributions in the proton has been 
observed in recent deep inelastic scattering^, FJ] and Drell-Yan experiments j|, || ||. 

There have been many theoretical attempts 0, || || 10, |ll]| trying to find the ori- 



gins for this asymmetry. It is believed 0, that the asymmetry cannot be produced 
from perturbative QCD and mesonic degrees of freedom play an important role for 
the effect. 
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In this paper we follow a new idea|12|] for the origin of the light flavor sea quark 
asymmetry to reproduce the recent observed d(x) — u(x) distribution @ with a sim- 
ple statistical model. The basic idea in Ref. 0] is rather simple: while sea quark- 
antiquark pairs are produced flavor blindly by gluon splitting, u quarks have larger 
probability to annihilate than d quarks due to the fact that there are more u quarks 
than d quarks in the proton, which hence causes the asymmetry. Taking proton as an 
ensemble of quark-gluon Fock states [0, [14]] and using the principle of detailed bal- 
ance for transitions between various Fock states through creation or annihilation of 
partons, the probabilities pij± of finding the quark-gluon Fock states \{uud}{i, j, k}) 
have been obtained and given in Table 2 of Ref. [12|, with i,j, k the number of uu 
pairs, the number of dd pairs, the number of gluons, respectively. With the density 
matrix Pij 7 k for the quark-gluon Fock states \{uud}{i, j, k}) , the sea-quark flavor 
asymmetry was calculated as d— u ph 0.124, which is in surprisingly agreement with 
the experimental data d — u = 0.118 ± 0.012. Encouraged by this success, here we 
want to extend the model to calculate the Bjorken-x distribution of partons to study 
the x-dependence of the flavor asymmetry in the nucleon sea, d(x) — u(x), which 
has recently been well measured by the FNAL E866/NuSea Collaboration^], [|. 

For a quark-gluon Fock state \{uud}{i, j, k}) , the total number of partons is 
n — 3 + 2i + 2j + k. If the n partons were free particles without mutual interactions, 
their momentum distribution dp^{pi, ■ ■ ■ ,p n ) would simply follow the n-body phase 



space, i.e. 



■ ■ 



n n j3 

d P F n {vu ■■■,Pn) = d<S> n (P; Pl , ■■-, Pn ) = 5\P - J>) j^ Wi (1) 
with P the 4-momentum of the proton, pi, ■ ■ ■ ,p n the 4-momenta of n partons. If 



we ignore the mass of partons, then E,i = \Jp 2 + mf = \pl\ = Pi and we have 

n n P./7P-WO- n n T7.rlI7.rlO- 

*f (p., ■■•.?»)= s'(p -gi*) n = s i (p - ep.) n (2) 

However, we know that partons are not free particles and are confined in the proton. 
In potential picture, partons are almost free only near the center of the proton and 
their momenta decrease when moving away from the center. Partons with smaller 
momenta at the center stay longer close to the center while partons with larger 
momenta at the center stay shorter around the center. So compared with total free 
particles, partons in the proton should have larger probability for smaller momenta. 
In this paper we assume the n-parton momentum distribution dp n (pi, ■ ■ ■ ,p n ) to be 

dpn{Pl,- ■ ■ ,Pn) = T-rn 1 p d^ n {P;p 1 ,---,p n ) 

= 8\P- Y /Pl )T]^^ = S\P- v Pi m^=. (3) 
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This is equivalent to assuming equal probability for any energy configuration (E\, ■ ■ ■ , E n ) 
of n partons in the proton. 

With the n-parton momentum distribution dp n (pi, ■ ■ ■ ,p n ) and taking the light 
front formula for the Bjorken-x 

-^parton Pz parton 



X 



M, 



(4) 



proton 



we can easily get the x-distribution of partons, p n (x), for an n-parton Fock state of 
the proton by using a simple Monte Carlo simulation with a Monte Carlo event gen- 
erator program called GENEV from the CERN computer program library (CERN- 
LIB). The p n (x) is normalized as Jq dxp n (x) = 1. 

For a n-parton Fock state \{uud}{i, j, k}), the total number of partons is n 
including u, d, u, d and g with parton number of 2 +i, 1 + j, i, j and k, respectively. 
Then the ^-distribution for each kind of partons, n^^x), dij,k( x )-> ny^x), dij,k( x ) 
and gij,k{x), are 



dij t k \X) 
9i,j,k {.•!') 



p n {x)(2 + i), 

p n (x)(l+j), 
p n (x)i, 
Pn(x)j, 
p n (x)k. 



(5) 
(6) 
(7) 
(8) 
(9) 



Summing up parton x-distribution of all possible Fock states \{uud}{i, j, k}) ac- 
cording to their weights (pij k) i n Table 2 of Ref . [pi to get the parton x-distribution 
of proton, we have 



U ( X ) = Yl Pi,j,k U i,j,k( X )i 

u( x ) = Y Pi,j,kUi,j,k(x), 
d( x ) = Y Pi,j,kdi,jA X )> 

9( x ) = Y P'.j k!l, .iA x )- 

i,3,k 



which satisfy normalization condition 



x 



u(x) + d(x) + u(x) + d(x) + g(x) 



dx 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 



Our u(x) and d(x) include both valence and intrinsic sea quarks which are identical 
and not distinguishable in our approach. The x-distribution of valance quarks (u v , 
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d v ) can be easily obtained as 



u v (x) 
d v (x) 



u(x) — u(x), 
d(x) — d(x). 



In addition, we have 



n 



E 



1 r 



u(x) + d(x) + u(x) + d(x) + g(x) dx = 5.57, 
M Proton 0.938GeV 



n 



5.57 



0.168GeV, 



(16) 
(17) 



(18) 
(19) 



where n is the average number of partons in proton and E is the average energy of 
partons in proton. 
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Figure 1: The parton density / (/ = u v , d v , u, d, g) of our model at a scale 
//o ~ 0.168 GeV. The f(x) in our model is simulated using Monte Carlo with- 
out any parameter. All we need in our model are principle of detailed balance and 
assumption of equal probability for every energy configuration of an n-parton Fock 
state. 

From above equations, we get the parton x-distribution densities f(x, (j,q) (/ = 
u, d, u, d, g) as shown in Fig.p] with the scale /j,q ~ E = 0.168 GeV. 

The corresponding momentum distribution, x/(x,//q), is shown in Fig.|2|. As a 
comparison, the xf(x, Q 2 ) distribution of GRVHH at some higher Q 2 is shown in 
FigfJ. It will be interesting to check whether with some QCD evolution equation 
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Figure 3: For a comparison, the xf(x) of GRV|]T3[ at higher Q 2 . Left: input densities 



xf (f = u v , d v , u, d, g) at Q 2 



0.26 GeV 2 and Q 2 



^NLO 



0.40 GeV 2 , at 



which scales strange sea s = s vanishes; Right: the evolved results at Q 2 = 5 GeV 2 
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our xf(x,Ho) distribution could evolve into a momentum distribution xf(x,Q 2 ) at 
higher Q 2 to be similar to the GRV's. 

The quarks and gluons in the Fock states are the "intrinsic" partons of the 
proton, multi-connected non-perturbatively to the valence quarksjnj. Such partons 
are different from "extrinsic" partons generated from the QCD hard bremsstrahlung 
and gluon-splitting as part of the lepton scattering interaction. Partons measured 
at certain Q 2 by experiments include both "intrinsic" and "extrinsic" ones. Since 
"extrinsic" partons are generated flavor blindly, the light flavor sea quark asymmetry 
is mainly due to "intrinsic" partons and is practically Q 2 independent, although 
correlations between "intrinsic" and "extrinsic" partons can cause some small Q 2 
dependence for the asymmetry. Experimental data at various Q 2 values also show 
little Q 2 dependence]]], ||. Hence we can compare our model prediction of d(x) —u(x) 
at the low scale fJLo directly with experimental data at higher scales, as shown in 
Fig.^. Our prediction is in good agreement with recent experiment data of FNAL 
E866/NuSea@] at Q 2 = 54 GeV 2 and HERMES |] at < Q 2 >= 2.3 GeV 2 . 




Figure 4: Comparison of measured d(x) — u(x) to prediction of our model at scale 
Ha « 0.168GeV. The FNAL E866/NuSea results, scaled to fixed Q 2 = 54GeV 2 , are 
shown as the circles; HERMES results of < Q 2 >= 2.3GeV 2 are shown as squares. 
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In summary, following Ref. |12j we take proton as an ensemble of quark-gluon 
Fock state, i.e., \p) = J2i,j,k c i,j,k\i uu d}{hj,k}), with pij k = | Ci,j,fc | 2 determined 
by the principle of detailed balance. By further assuming equal probability for 
any energy configuration (E±, ■ • • , E n ) of n-parton Fock state in the proton, we get 
parton x-distribution functions at a scale of /io ~ E = 0.168 GeV with a Monte 
Carlo simulation of a simple statistical model. The corresponding light flavor sea 
quark asymmetry d(x) — u(x) reproduces the recent experiment data quite well. This 
is a further support of the new origin of the light flavor sea quark asymmetry as a 
pure statistical effect due to the fact that there are more u-quarks than (i-quarks in 
the proton. 
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